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Abstract 

The quark-meson coupling model, based on a mean field description of non- 
overlapping nucleon bags bound by the self-consistent exchange of a, u and p 
mesons, is extended to investigate the properties of finite nuclei. Using the Born- 
Oppenheimer approximation to describe the interacting quark-meson system, we 
derive the effective equation of motion for the nucleon, as well as the self-consistent 
equations for the meson mean fields. The model is first applied to nuclear matter, 
after which we show some initial results for finite nuclei. 
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1 Introductory remarks 

The nuclear many-body problem has been the object of enormous theoretical attention 
for decades. Apart from the non-relativistic treatments based upon realistic two-body 
forces jlj, ||, there are also studies of three-body effects and higher ||. The importance of 
relativity has been recognised in a host of treatments under the general heading of Dirac- 
Brueckner 0, |5|, [J. This approach has also had considerable success in the treatment of 
scattering processes @. At the same time, the simplicity of Quantum Hadrodynamics 
(QHD) H [J has led to its widespread application, and to attempts to incorporate the 
density dependence of the couplings [TO, ITTJ that seems to be required empirically. 



One of the fundamental, unanswered questions in this field concerns the role of sub- 
nucleonic degrees of freedom (quarks and gluons) in determining the equation of state. 
There is little doubt that, at sufficiently high density (perhaps 5 — 10po, with p the 
saturation density of symmetric nuclear matter), quarks and gluons must be the correct 
degrees of freedom and major experimental programs in relativistic heavy ion physics are 



either planned or underway [12], [13| to look for this transition. In order to calculate the 



properties of neutron stars |14| one needs an equation of state from very low density to 



many times po a t the centre. A truly consistent theory describing the transition from 
meson and baryon degrees of freedom to quarks and gluons might be expected to incor- 
porate the internal quark and gluon degrees of freedom of the particles themselves. Our 
investigation may be viewed as a first step in this direction. We shall work with the quark- 



meson coupling (QMC) model originally proposed by one of us |T_5| and since developed 



extensively |T6| , |TT1 . Related work has been carried out by a number of groups [JJ| . 



Within the QMC model [|l5j the properties of nuclear matter are determined by the 
self-consistent coupling of scalar (a) and vector (u) fields to the quarks within the nucleons, 
rather than to the nucleons themselves. As a result of the scalar coupling the internal 
structure of the nucleon is modified with respect to the free case. In particular, the small 
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mass of the quark means that the lower component of its wave function responds rapidly 
to the a field, with a consequent decrease in the scalar density. As the scalar density 
is itself the source of the a field this provides a mechanism for the saturation of nuclear 
matter where the quark structure plays a vital role. 

In a simple model where nuclear matter was considered as a collection of static, non- 
overlapping bags it was shown that a satisfactory description of the bulk properties of 
nuclear matter can be obtained []15 , 16| . Of particular interest is the fact that the extra 
degrees of freedom, corresponding to the internal structure of the nucleon, result in a 
lower value of the incompressibility of nuclear matter than obtained in approaches based 
on point-like nucleons - such as QHD ||. In fact, the prediction is in agreement with the 
experimental value once the binding energy and saturation density are fixed. Improve- 
ments to the model, including the addition of Fermi motion, have not altered the dominant 
saturation mechanism. Furthermore, it is possible to give a clear understanding of the 
relationship between this model and QHD Jl6j and to study variations of hadron proper- 



ties in nuclear matter ||17|| . Surprisingly the model seems to provide a semi- quantitative 
explanation of the Okamoto-Nolen-Schiffer anomaly JT^| when quark mass differences are 
included |[20|| . A further application of the model, including quark mass differences, has 
suggested a previously unknown correction to the extraction of the matrix element, V u d, 
from super-allowed Fermi beta-decay ||21|| . Finally the model has been applied to the case 
where quark degrees of freedom are undisputedly involved - namely the nuclear EMC 
effect |3 . 

Because the model has so far been constructed for infinite nuclear matter its applica- 
tion has been limited to situations which either involve bulk properties or where the local 
density approximation has some validity. Our aim here is to overcome this limitation 
by extending the model to finite nuclei. In general this is a very complicated problem 
and our approach will be essentially classical. Our starting point will be exactly as for 
nuclear matter. That is, we assume that on average the quark bags do not overlap and 
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that the quarks are coupled locally to average a and 00 fields. The latter will now vary 
with position, while remaining time independent - as they are mean fields. For deformed 
or polarised nuclei one should also consider the space components of the vector field. In 
order to simplify the present discussion we restrict ourselves to spherical, spin-saturated 
nuclei. 

Our approach to the problem will be within the framework of the Born-Oppenheimer 
approximation. Since the quarks typically move much faster than the nucleons we assume 
that they always have time to adjust their motion so that they are in the lowest energy 
state. In order to account for minimal relativistic effects it is convenient to work in 
the instantaneous rest frame (IRF) of the nucleon. Implicitly one then knows both the 
position and the momentum of the nucleon, so that the treatment of the motion of the 
nucleon is classical - at least, as long as the quarks are being considered explicitly. The 
quantisation of the motion of the nucleon is carried out after the quark degrees of freedom 
have been eliminated. 

In Sect. H] we show how the Born-Oppenheimer approximation can be used to treat the 
quark degrees of freedom in a finite nucleus. In Sect. ^ we derive the classical equation of 
motion for a bag in the meson fields, including the spin-orbit interaction which is treated 
in first order in the velocity. We then quantize the nucleon motion in a non relativistic 
way. The self consistent equations for the meson fields are derived in Sect. [|. In Sect. [5] 
we summarize the model in the form of a self-consistent procedure. The use of this 
non relativistic formulation is postponed to future work. To allow a clear comparison 
with QHD we propose a relativistic formulation in Sect. |^, which is then applied to 
nuclear matter in Sect. [?]. Some initial results for finite nuclei are also presented. Sect. § 
summarises our main results and identifies directions for future work. 
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2 The Born-Oppenheimer approximation 

In what follows the coordinates in the rest frame of the nucleus (NRF) will be denoted 
without primes: (t,r). In this frame the nucleon follows a classical trajectory, R(t). 
Denoting the instantaneous velocity of the nucleon as v = dR/dt, we can define an 
instantaneous rest frame for the nucleon at each time t. The coordinates in this IRF are 

(t>,r>y- 

r L = r^cosh£ + £'sinh£, 

f± = rV, (1) 

t = t'cosh£ + r^sinh^, 

where tl and r± are the components respectively parallel and transverse to the velocity 
and £ is the rapidity defined by tanh£ = \v(t)\. 

Our assumption that the quarks have time to adjust to the local fields in which the 
nucleon is moving is exact if the fields are constant - i.e. if the motion of the nucleon has no 
acceleration. It is, of course, very important to examine the validity of the approximation 
for a typical nuclear environment. For this purpose we take the nucleon motion to be 

— * 

non-relativistic. Assume that at time the nucleon is at -Ro- After time t, assuming t 
small enough, we have 

R(t) =R + v t + h 2 a , (2) 

where a (= F/M N = —VV/M N ) is the acceleration and M N is the free nucleon mass. 
We shall take the potential to be a typical Woods-Saxon form with depth V ~ —50 MeV, 
surface thickness a ~ 2 fm and radius Ra ~ 1.2 A 1 ^ 3 fm. The maximum acceleration 
occurs at it! = Ra and takes the value 

R 

Umax = ^0—7^—- (3) 



Therefore, in the IRF we have 



^W-S + i^A,, (4) 



and, in the worst case, relative to the size of the nucleon itself {Rb), the departure from 
a fixed position, is 



fl'(t) - R' 



„ *, (5 ) 

n n 2 aM N R B 80 v ' 

with t in fm. Thus, as long as the time taken for the quark motion to change is less 

than ~ 9 fm, the nucleon position in the IRF can be considered as unchanged. Since the 

typical time for an adjustment in the motion of the quarks is given by the inverse of the 

typical excitation energy, which is of order 0.5 fm, this seems quite safe. 

As we have just seen, it is reasonable to describe the internal structure of the nucleon 

in the IRF. In this frame we shall adopt the static spherical cavity approximation to the 



MIT bag for which the Lagrangian density is 

Co = tf{i-fdp - m q )q' - BV, for \u'\ < R B , (6) 

with B the bag constant, Rb the radius of the bag, m q the quark mass and u' the position 
of the quark from the center of the bag (in the IRF): u ' = f ' — R' . We shall denote as u' 
the 4-vector (t',u'). The field q'(u' ,u') is the quark field in the IRF, which must satisfy 
the boundary condition 

(1 + ij ■ u')q' = 0, sX\u'\=Rb- (7) 

Next we must incorporate the interaction of the quarks with the scalar (a) and vector 
(u) mean fields generated by the other nucleons. In the nuclear rest frame they are self- 
consistently generated functions of position - a(r) and u(r). Using the scalar and vector 
character of these fields, we know that in the IRF their values are 

criRF(t',u') = cr(f), 

uj IRF (t',u') = ^(f)cosh£, ( 8 ) 
£uiRF(t' ,u') = — uj(f)vsmh£, 

and the interaction term is 

C Z = glq'q\u')a IRF (u') - g q J^q'(u')u? RF (u'), (9) 



where g q a and g% are the quark- meson coupling constants for a and uj, respectively. Apart 
from isospin considerations, the effect of the p meson can be deduced from the effect of 
the uj. Thus we postpone its introduction to the end of this section. 

Since we wish to solve for the structure of the nucleon in the IRF we need the Hamil- 
tonian and its degrees of freedom in this frame. This means that the interaction term 
should be evaluated at equal time t' for all points u' in the bag. Suppose that at time t' 
the bag is located at R! in the IRF. Then, in the NRF it will be located at R at time T 
defined by 

R L = i^coshf + i , sinhf, 

R± = R' ± , (10) 
T = t'cosh^ + i^sinh^. 

For an arbitrary point r ' (= u' + R!) in the bag, at the same time t', we have an analogous 
relation 

r L = r^cosh£ + t'sinh£, 

rj_ = rV, (11) 
t = t'cosh£ + r^sinh£, 

from which we deduce 

r L = R' L cosh ^ + t' sinh^ + u' L cosh ^, 

= R L + M' L cosh^, (12) 
f± = R± + u± ' . 

This leads to the following expression for the a field 

<JiR F (t',u') =a(R L (T)+u' L cosh^R ± (T)+u ± '), (13) 

with a corresponding equation for the uj field. 



The spirit of the Born-Oppenheimer approximation is to solve the equation of motion 
for the quarks with the position R(T) regarded as a fixed parameter. In order to test 
the reliability of this approximation we consider a non-relativistic system and neglect 
finite-size effects. That is, we take 

a IRF (t',u')~a(R(T)). (14) 

As we noted earlier, the typical time scale for a change in the motion of the quark is t ~ 
0.5 fm. During this time the relative change of a due to the motion of the bag is 

^ = v-R a -r. (15) 
a a 

It is reasonable to assume that a roughly follows the nuclear density, and as long as this is 
constant, Aa vanishes and the approximation should be good. The variation of the density 
occurs mainly in the surface where it drops to zero from p (the normal nuclear density) 
over a distance d of about 2 fm. Therefore we can estimate \o' ja\ as approximately 1/d 
in the region where p varies. The factor v ■ R depends on the actual trajectory, but as 
a rough estimate we take v ■ R ~ 1/3. That is, we suppose that the probability for v 
is isotropic. For the magnitude of the velocity we take kpjM^ with k F = 1.7 fm -1 the 
Fermi momentum. With these estimates we find 

Aa 0.36-0.5 

~ o o 3 %> ( 16 ) 

a 6-1 

which is certainly small enough to justify the use of the Born-Oppenheimer approximation. 
Clearly this amounts to neglecting terms of order v in the argument of a and uj. In order 
to be consistent we therefore also neglect terms of order v 2 - that is, we replace u' L cosh £ 
by u' L . 
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3 Equation of motion for a bag in the nuclear field 
3.1 Leading term in the Hamiltonian 

Following the considerations of the previous section, in the IRF the interaction Lagrangian 
density takes the simple, approximate form: 

Ci = g g Jq\t>, u')a{R + u') - g q J[ lo cosh £ + 7 ■ v sinh £]g'(t', u')u(R + u'). (17) 

The corresponding Hamiltonian is 

H = / du'q'[-i^- V + m q - g q a a(R + u') 
Jo 

+g^uj(R + u')^° cosh^ + ^ ■ vsmh^)]q'(t',u') + BV, (18) 
while the momentum is simply 

P= / du'q'^-iVW. (19) 
Jo 

As the a and u fields only vary appreciably near the nuclear surface, where R ^> \u'\ 
(since \u'\ is bounded by the bag radius), it makes sense to split H into two parts: 

H = H + H 1 , (20) 
rR B 

H = / du'q'[-tf- V + m q -gla{R) 
Jo 

+g q u(R)( 1 ° cosh£ + 7 ■ v smh£)}q'(t' ,u') + BV, (21) 
rRB 

H x = / du'q'[-gl{a{R + u')-a{R)) 
Jo 

+gl(u(R + u')- ou(R))(7 cosh£ + 7 • v smh^W, u'), (22) 

and to consider Hi perturbation. 

Suppose we denote as 4> a the complete and orthogonal set of eigenfunctions defined 

by 

h(j) a (u') = (-i 7 °7- V + my)(j) a (u'), 



■"■Only the quark degrees of freedom are active. The nucleon position and velocity are parameters as 
discussed earlier. 
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;i + *7-« )0 Q (« ) 

/ du'^<f = S a/3 , 
Jo 



Rb 

0, at \u'\ = Rb, 



(24) 
(25) 



with {a} a collective symbol for the quantum numbers and m* a parameter. Here we 
recall the expression for the lowest positive energy mode, Om , with m the spin label: 

j (xu'/R B ) 



where 



Xr, 



i(5 q a ■ u' 'ji{xu' I 'R B ) ) VlP 



(26) 



AT 2 



x 2 + {m*R B ) 2 , /3 q 



tt - m*R B 



2R B j 2 (x)[Q (Q - 1) + m*i? B /2]/x^ 



For this mode, the boundary condition at the surface amounts to 



jo(x) = P g ji(x). 



We expand the quark field in the following way 



q'{t',u') = Y,e-^ U 'ct ) a {u')b a {t') 1 



(27) 
(28) 



(29) 



(30) 



with k chosen as 

k = g q lo (R)v smb. (31) 

in order to guarantee the correct rest frame momentum for a particle in a vector field. 
Substituting into the equation for H we find 



a 



JIB a 

a ap 



+N q g q w u{R) coshf + BV 
P = £>| - *V|/9)6t6^ - N q k, 



(32) 
(33) 
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with the notation 

i-Rb 



A\p) = du' <j> a Hu')Af{u') and N q = Yb)Ja a . (34) 
Jo „ 



Choosing m* = m q — g^a(R) (in which case the frequency Q a and the wave function 
4> a become dependent on R through a) we get for the leading part of the energy and 
momentum operator 

H j RF = ^^) 6 t 6a + W + i v^ M ^ )coshe) (35) 



a 



Rb 



P IRF = Y( a \-W\f3)b%-N q g^(R)vsmhZ. (36) 

a/3 

If we quantize the b a in the usual way, we find that the unperturbed part of H is di- 
agonalised by states of the form \N a ,Np, • • •) with A^ a the eigenvalues of the number 
operator b^ a b a for the mode {&}■ According our working hypothesis, the nucleon should 
be described by three quarks in the lowest mode (a = 0) and should remain in that 
configuration as R changes. As a consequence, in the expression for the momentum, the 
contribution of the gradient acting on <fi averages to zero by parity and we find that the 
leading terms in the expression for the energy and momentum of the nucleon in the IRF 



arc: 



E T RF = M* N (R)+3g«oj(R) C osht (37) 
P IRF = -3gluj(R)v sinh£, (38) 



with 



MUR) = ^^ + BV. (39) 

Since we are going to treat the corrections to leading order in the velocity, they will 
not be affected by the boost back to the NRF. Therefore we can apply the Lorentz 
transformation to Eqs. ( |37|) and (|38| ) to get the leading terms in the energy and momentum 
in the NRF. Since there is no possibility of confusion now, we write the NRF variables 
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without the NRF index. The result is 

E = M* N (R) cosh £ + 3gZu{R), (40) 
P = M* N {R)v sinh£, (41) 

which implies 



E = \]Mf{R) + P 2 + 3g*u(R). (42) 

— * 

At this point we recall that the effective mass of the nucleon, M^(R), defined by 
Eq.(p9|) does not take into account the fact that the center of mass of the quarks does 
not coincide with the center of the bag. By requiring that all of the quarks remain in the 
same orbit one forces this to be realized in expectation value. However, one knows that 
the virtual fluctuations to higher orbits would decrease the energy. This cm. correction is 
studied in detail in the Appendix, where it is shown that it is only very weakly dependent 
on the external field strength for the densities of interest. For the zero point energy due 
to the fluctuations of the gluon field, we assume that it is the same as in free space. Thus 
we parametrize the sum of the cm. and gluon fluctuation corrections in the familiar form, 
—zq/Rb, where Zq is independent of the density. Then the effective mass of the nucleon 
in the nucleus takes the form 

M* N (R) = 3n ° { ^- Z ° +BV, (43) 
Kb 



and we assume that the equilibrium condition is 

dM* N {R) 



0, (44) 



dRs 

which is the usual non-linear boundary condition. This is again justified by the Born- 
Oppenheimer approximation, according to which the internal structure of the nucleon has 
enough time to adjust to the varying external field so as to stay in its ground state. 

The parameters B and Zq are fixed by the free nucleon mass (M/y = 939 MeV) using 
Eqs.(|43|) and ([f4]) applied to the free case. In the following we keep the free bag radius, 
R B , as a free parameter. The results are shown in Table [|. 
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Table 1: SV*(MeV) and z for some bag radii using m q = 5 MeV. 





0.6 


0.8 


1.0 


£1/4 


210.9 


170.0 


143.8 


^0 


4.003 


3.295 


2.587 



3.2 Corrections due to H\ and Thomas precession 

We estimate the effect of Hi in perturbation theory by expanding a(R+u') and u(R + u') 
in powers of u' and computing the effect to first order. To this order several terms give 
zero because of parity and one is left with 

((Ofl^KO) 3 ) =^E((0) 3 |&l(0) 3 )(«|7 7-^'sinhe|/3) ■ V R u(R). (45) 

a/3 

Here we need to be more precise about the meaning of the labels a, (3. Suppose we set 
{a} = {0,m a } with m a the spin projection of the quark in mode {0}, then we find 

((0)m a |7°7 • vu'smh £1(0)771/?) = —I(R)(m a \^-\mp) x f)sinh£, (46) 



with 



2 



4 rRB 

I(R) = -Af 2 / du'u' 3 Jo (xu'/R B )(3 q ji(xu'/R B ), (47) 



3 Jo 
R B ( 4fi + 2m*R B - 3 



(48) 



3 \2Q (Q - 1) + m*R B J 

— * — * 

The integral, I(R), depends on R through the implicit dependence of R B and x on 
the local scalar field. Its value in the free case, Iq, can be related to the nucleon isoscalar 
magnetic moment: J = 3fi s /M N with /x s = /i p + // n and \i p = 2.79, // n = —1.91 the 
experimental values. By combining Eqs.fl45D, (TO) and (fH"D we then find 



((o) 3 !^) 3 ) = /i. ^ m~nR ^ x dshlH (49) 

^ s I M N M%{R)R \dR ^ 1 7 1 ; 
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with S the nucleon spin operator and L its angular momentum. 

This spin-orbit interaction is nothing but the interaction between the magnetic mo- 
ment of the nucleon with the "magnetic" field of the u seen from the rest frame of the 
nucleon. This induces a rotation of the spin as a function of time. However, even if 
li s were equal to zero, the spin would nevertheless rotate because of Thomas precession, 
which is a relativistic effect independent of the structure. It can be understood as follows. 

Let us assume that at time t, the spin vector is S(t) in the IRF(t). Then we expect that, 
at time t + dt the spin has the same direction if it is viewed from the frame obtained by 
boosting the IRF(t) by dv so as to get the right velocity v(t + dt). That is, the spin looks at 
rest in the frame obtained by first boosting the NRF to v(t) and then boosting by dv. This 
product of Lorentz transformation amounts to a boost to v(t + dt) times a rotation. So, 
viewed from the lKF(t + dt), the spin appears to rotate. In order that our Hamiltonian be 
correct it should contain a piece H prec which produces this rotation through the Hamilton 



equations of motion. A detailed derivation can be found in Refs. p4| , 25| and the result is 



Hprec = X ^ • S. (51) 

The acceleration is obtained from the Hamilton equations of motion applied to the leading 
order Hamiltonian, Eq . (ff2|) . To lowest order in the velocity one finds 

f t -^ lM '^ ) + 3M (52) 

If we put this result into Eq. (|5"ID and add the result to Eq. (|50D , we get the total spin orbit 
interaction (to first order in the velocity) 

H prec + H l = V s . {R)S-L, (53) 

where 

V S . .(R) : 



d M* N (R)) + (1 - 2^.(3)) ( 4eWZm>{R) 



2M* N 2 (R)R[\dR JVV 7 v "\dR 

and 



(54) 



n I(R)M* N (R) 

1qM n 
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3.3 Total Hamiltonian for a bag in the meson mean fields 

To complete the derivation we now introduce the effect of the neutral p meson. The 
interaction term that we must add to Ci (see Eq. ([)])) is 

C P i = -g^l^iu')^ IRF {u'l (56) 

where p^ IRF is the p meson field with isospin component a and r a are the Pauli matrices 
acting on the quarks. In the mean field approximation only a = 3 contributes. If we 
denote by b(R) the mean value of the time component of the field in the NRF, we can 
transpose our results for the uj field. The only difference comes from trivial isospin factors 
which amount to the substitutions 

~> 9p^, Pv = Pp~ Pn, (57) 

where (with eigenvalues ±1/2) is the nucleon isospin operator. 

This leads to our final result for the NRF energy-momentum of the nucleon moving 
in the mean fields, a(R),u(R) and b(R): 

E = M$f(R) cosh £ + V(R), (58) 
P = M*(i?){) sinh£, (59) 



with 



V(R) = V C (R)+V S . .(R)S-L, (60) 



V C (R) = 3glu(R)+g^b(R), (61) 

1 - - T? 

V s .o.{R) = -^—[A a + (l-2p s r ]s (R))A UJ + (l-2p v rj v {R))^-A p \, (62) 



2M$(R)R L " " v ru '" x " 2 

K = ^M**(R), \> = Jr^MR), A p = ^9 q P KR). (63) 



Up to terms of higher order in the velocity, this result for the spin-orbit interaction 
agrees with that of Achtzehnter and Wilets |f26| who used an approach based on the non- 



topological soliton model. As pointed out in ||26|| , for a point-like Dirac particle one has 
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/jl s — 1 while the physical value is /i s = .88. Thus, in so far as the omega contribution 
to the spin-orbit force is concerned, the point-like result is almost correct. This is clearly 
not the case for the rho contribution since we still have /i v — 1 for the point-like particle 
while experimentally /i v = 4.7. 

3.4 Quantization of the nuclear Hamiltonian 

Until now the motion of the nucleon has been considered as classical, but we must now 
quantize the model. We do that here in the non-relativistic framework. That is, we 
consider a theory where the particle number is conserved and we keep only terms up to 
those quadratic in the velocity. (Here we drop the spin dependent correction since it 
already involves the velocity. It is reinserted at the end.) 

The simplest way to proceed is to realize that the equations of motion 

E = M^(R) cosh £ + V c (R), (64) 
P = M* N (R)v sinh£, (65) 

can be derived from the following Lagrangian 

L(R, v) = -M^(R)Vl-v 2 - V C (R). (66) 

Thus the expansion is simply 

L nr (R, v) = \m* n {R)v 2 - M* N (R) - V C (R) + <D(v*) (67) 

Neglecting the 0(v 4 ) terms we then go back to the (still classical) Hamilton variables and 
find (note that the momentum is not the same as before the expansion) 

H nr (R, P) = ——^ + Mtf(R) + V C (R). (68) 
2M N (K) 

If we quantize by the substitution P — > —iVr we meet the problem of ordering ambi- 
guities since now M%{R) and P no longer commute. The classical kinetic energy term 
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P 2 /2M^(R) allows 3 quantum orderings 

T x = P ■ AP, T 2 = AP 2 , T 3 = P 2 A, (69) 

where A = 1/(2M^ (R)). However, since we want the Hamiltonian to be hermitian the 
only possible combination is 

T = zT 1 + ^(T 2 + T 3 ), (70) 

with z an arbitrary real number. Using the commutation relations one can write 

T 2 + T 3 = 2T 1 -i^(V^i), (71) 

so that the kinetic energy operator 

T = T 1 -^(V%A), (72) 

is only ambiguous through the term containing the second derivative of A = 1/(2M^(R)). 
Since we have expanded the meson fields only to first order in the derivative to get the 
classical Hamiltonian it is consistent to neglect such terms. So our quantum Hamiltonian 
takes the form 

H nr (R, P)=P. —^-P + M* N (R) + V(R), (73) 
ZM N {K) 

where we have reintroduced the spin-orbit interaction in the potential V(R). The nuclear 
quantum Hamiltonian appropriate to a mean field calculation is then 

H nr = ]T H nr (R u P t ), P t = -iVi. (74) 

i=l,A 

The problem is manifestly self consistent because the meson mean fields, upon which H nr 
depends through and V, themselves depend on the eigenstates of H nr . 

4 Equations for the meson fields 

The equation of motion for the meson-field operators (<r, oj u , p u,a ) are 

d^a + mla = g«qq, (75) 
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dpW + mfa 1 ' = gfafq, (76) 
d^p^ + m 2 p p^ = fltf^q. (77) 

The mean fields are defined as the expectation values in the ground state of the nucleus, 
\A): 

(A\a(t,f)\A) = a(r), (78) 
(A\^(t,r)\A) = 6(y,0)u(f), (79) 
(A\p^(t,r)\A) = 6(y,0)6(a,3)b(f). (80) 

The equations which determine them are the expectation values of Eqs.([75j), Q7T3D and 
([77|). First we need the expectation values of the sources 

(A\qq(t,r)\A),(A\qYq(t,f)\A) and (A\qj v ^-q(t, r)\A) . (81) 

As before, we shall simplify the presentation by not treating the p meson explicitly until 
the end. In the mean field approximation the sources are the sums of the sources created 
by each nucleon - the latter interacting with the meson fields. Thus we write 

qq(t,r) = £ (qq(t,r)) u (82) 
i=l,A 

qi v q{t^ = E (<tfq(t,r))i, (83) 

i=l,A 

where (• • •)» denotes the matrix element in the nucleon i located at Ri at time t. According 
to the Born-Oppenheimer approximation, the nucleon structure is described, in its own 
IRF, by 3 quarks in the lowest mode. Therefore, in the IRF of the nucleon i, we have 

(3W,0>< = Sj2t l ' m (u')^' m (u') = Ss t (u'), (84) 

m 

(?' 7 Y(t',r")>* = 35(u,0)Y,<p]°' m (u , )^ m (u')=35(u,0)w l (u'), (85) 

m 

where the space components of the vector current gives zero because of parity. At the 
common time t in the NRF we thus have 

R' i L = Ri tL cosh& - t sinh^, R' i ± = R i>± , (86) 
r' L = tl cosh£, — t sinh^j, r±' = r±. (87) 
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Therefore 

u i,L = ( r L ~ Ri,l) cosh&, u i:± ' = f± - R i:± , 
and from the Lorentz transformation properties of the fields we get 

— * 

(qq(t,r))i = 3s;((r L - R ijL ) cosh&, r± - R it ±), (89) 

{q r y°q{t,r)) i = 3w;((r L - R ijL ) cosh r ± - R ii± ) cosh (90) 

(q^q(t,r))i = 3iVi((r L - R^l) cosh&, rj_ - -R^^sinh^. (91) 

These equations can be re-written in the form 

(??(f,r))i = (^3( cosh ^) _1 Jdke^-^ S(k,Ri), (92) 
(57°?(*,*0>i = j^Jdke^-^W&Ri), (93) 
(fffa(*,*0>i = j^Vijdke^-^Wik,^), (94) 



with 



5ft = / due- l( ^^ +kLUL/cosh ^ Si (u), (95) 

W(jfe,flj) = J due~ i{% ^ +kLULlcosh ^ w t {u). (96) 

Finally, the mean field expressions for the meson sources take the form 

(A\qq{t^\A) = ^l^e^l^^osh^^e-^^ft^)^), (97) 

(A\vfq(t,i>)\A) = j^fdke a *(A\^e^W(k,Ri)\A), (98) 

(A\q^q(t,r)\A) = 0, (99) 

where the last equation follows from the fact that the velocity vector averages to zero. 
To simplify further, we remark that a matrix element of the form 

(A\J2e- iU - Ri ---\A) (100) 



20 



is negligible unless k is less than, or of the order of, the reciprocal of the nuclear radius. 
But in Eqs. (^) and ( p6|) k is multiplied by u which is bounded by the nucleon radius. 
Hence, if we restrict the application of the model to large enough nuclei, we can neglect 
the argument of the exponential in Eqs. (pop and (||J). The evaluation of the correction 
to this approximation will be postponed to a future work. 

We now define the scalar, baryonic and isospin densities of the nucleons in the nucleus 

by 

Ps(r~) = {A^ ^fiA r-Rim, (101) 
p B (r) = (A| AM, (102) 

i 

P 3 (r) = (A\J2\Kr-Rd\A), (103) 

where we have used Eq. (|58"D to eliminate the factor (cosh£j) -1 . Note that the definition 
of the scalar density makes sense because, in mean field approximation, each nucleon is 
moving in an orbital with a given energy. The meson sources then take the simple form 

(A\qq(t,r)\A) = 3S(f)p s (r), (104) 
(A\qYq(t,r)\A) = 36(v,0)p B (r), (105) 
{A\qrf^q{t,¥)\A) = 5(u,0)5(a,3)p 3 (r), (106) 

where we have deduced the source of the p from that for the u. We have used the notation 

S{r) = S{6,f) = fdusriu), (107) 

(108) 



n /2 + rn*R B {n -l) 



n (n -l)+m*R B /2' 

where the subscript r reminds us that the function s?(u) must be evaluated in the scalar 
field existing at r. 

Since their sources are time independent and since they do not propagate, the mean 
meson fields are also time independent. So by combining Eqs.(|75|) to ( p0|) and Eqs.( |104| ), 
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( |105| ) and ( |106| ), we get the desired equations for a(r) } uj(r) and b(r): 

(-V 2 r + ml)a(r) = 9<J C(r)p s (r), (109) 
(-V r 2 + mXf) = 9ujPB {r), (110) 
(-V r 2 + mJ)6(f) = g p p 3 (f). (Ill) 

where the nucleon coupling constants and C are defined by 

</, = 3<#% = 0), g u = 3&, g P = g% C(r) = S(r)/S(a = 0). (112) 

For completeness we recall that the mean fields carry the following energy 

E meson = }_ J ^[(^2 + ^2 _ (y^2 _ ^2 _ (y^ _ ^ 

5 The problem of self-consistency 

Our quark model for finite nuclei is now complete. As the main equations are scattered 
through the text, it is useful to summarize the procedure which we have developed as 
follows: 

1. Choose the bare quark mass, m q , and adjust the bag parameters, B and Zq, to fit 
the free nucleon mass and its bag radius (see Eq . (^3|) and Table [I]). 

2. Assume that the coupling constants and the masses of the mesons are known (see 
Table | in Sec. ^D- 



3. Evaluate the nucleon properties, 1(a) (Eq.(pE8|)) and S(a) (Eq.( |108| )), for a range of 
values of a (see also Eqs. ( |134| ) and (|135| ) in Sec. ffTT] ). 



4. Guess an initial form for the densities, p s (^),Ps(^*) and P3(r), in Eqs.( |101| ), ( |102| ) 



and (JLOJ). 

5. For p s (r) fixed, solve Eq. (|109|) for the a field. 
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6. For ps(r) and p%{f) fixed, solve Eqs. flllOQ and ( |111| ) for the u and p fields. 



7. Evaluate the effective mass M^(f) and the potential V(r) according to Eqs. (|43|) and 
(|60D . The bag radius at each point in the nucleus is fixed by Eq.([H]). For practical 
purposes it is useful to note that M^(r) and C(r) depend only on the value of the 
a field at r and that a simple parametrization, linear in g a a, works extremely well 
at moderate densities (see Eqs.( |134j ) and (|136|) in Sec. |7~T|) . So one does not need 



to solve the bag equations at each point when solving the self-consistent nuclear 
problem. 

8. Solve the eigenvalue problem defined by the nuclear Hamiltonian Eq.(|74|) and gener- 
ate the shell model from which the densities p s (r), pb{t) and P3(r) can be computed 
according to Eqs. (TO, ( TO and (|103D - 

9. Go to 5 and iterate until self-consistency is achieved. 

This procedure has to be repeated for each nucleus, which certainly implies considerable 
numerical work. We plan to study the implications of this non-relativistic formulation in 
a future work. 



6 Relativist ic formulation 

Here we attempt to formulate the model as a relativistic field theory for the nucleon in 
order to have a direct comparison with the widely used QHD. We make no attempt to 
justify the formulation of a local relativistic field theory at a fundamental level because this 
is not possible for a composite nucleon. Our point is that, in the mean field approximation, 
QHD has had considerable phenomenological success. We therefore try to express our 
results in this framework. The idea is to write a relativistic Lagrangian and to check that, 
in some approximation, it is equivalent to our non-relativistic formulation. 
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As shown earlier, our basic result is that essentially the nucleon in the meson fields 
behaves as a point like particle of effective mass M^(a(r)) moving in a potential g UJ uj{f). 
To simplify we do not consider the p coupling and therefore we shall only apply the model 
to iV = Z nuclei. 

As already pointed out the spin-orbit force of the u is almost the same as the one of 
a point like Dirac particle. So a possible Lagrangian density for this system is 

C = ^7 • dip - M* N (a)ipil) - g^ipl^ip + C mesonS} (114) 

where ip, a and are respectively the nucleon, a and u field operators. The free meson 
Lagrangian density is 

Cmesons = \{d^a - m\a 2 ) - U a uo u {d^ u - + ^m^^. (115) 

The comparison with QHD is straightforward. If we define the field dependent coupling 
constant g a ((j) by 

M* N {&) = M N - g a {&)&, (116) 

it is easy to check that g a (a = 0) is equal to the coupling constant g a defined in Eq.( |112| ). 
Then we write the Lagrangian density as 

C = iijrf ■ dip - M N ipip + g a (a)aipip - g^ipj^ip + C mesonsi (117) 

and clearly the only difference from QHD lies in the fact that the internal structure of the 
nucleon has forced a known dependence of the scalar meson-nucleon coupling constant on 
the scalar field itself. Note that this dependence is not the same as the one adopted in the 
density dependent hadron field theory [jTTJ where the meson-nucleon vertices are assumed 
to depend directly on the baryonic densities. 

In the mean field approximation, the meson field operators are replaced by their time 
independent expectation values in the ground state of the nucleus: 

a(t,r) - a(f), u>»(t,r) -> 8(jji,0)u(r), (118) 
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and variation of the Lagrangian yields the Dirac equation 

(i 7 • d - M* N {a) - g u >y<w)ij> = 0, (119) 

as well as the equations for the meson mean fields 

(-V* + m 2 >(r) = - QUf* (a)\ (A\^\A) t (120) 
(-V 2 r + mlMr} = g u (A\tfif>tf\A). (121) 



Using the Foldy-Wouthuysen transformation (see Ref. [26] for the details), one can 
show that the Dirac equation ( |119| ) gives back our non-relativistic, quantum Hamiltonian 
(without the p coupling) under the following conditions: 

• only terms of second order in the velocity are kept, 

• second derivatives of the meson fields are ignored, 

• the fields are small with respect to the nucleon mass, 

• the difference between p s = .88 and p s {point) = 1 can be neglected. 

The fact that fi v = 4.7 is very different from p v (point) = 1 prevents one from using this 
simple scheme for the p coupling. The spin-orbit potential would be completely different 
from the one obtained in Eq.fl62"l). Since the p is not required in the treatment of symmetric 
nuclei we postpone to future work its introduction in the relativistic framework. 

In the same approximations one finds that (A\ipif)(r)\A) and (A\i/)'ip(r)\A) are respec- 
tively equal to the previously defined scalar source, p s (Eq. fllOlp ), and vector source, 
Pb (Eq. (|102|) ). (In nuclear matter the approximation is in fact exact). Therefore the 
equations (|120|) and ( |121|) for the meson fields reproduce our previous results given in 
Eqs. (|109p and ( |110|) provided the relation 



C(a)g a (a = 0) = -^-M^(a), (122) 
= ^(^(a)a), (123) 
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is satisfied, which can be checked explicitly from the definition of Mjy-(er) and C(a). 

Within the limits specified above, the Lagrangian density given in Eq.( |114j) looks 
acceptable and we shall proceed to study its numerical consequences. 



7 Applications 

Before turning to finite nuclei, we first explain briefly how the general formalism applies 
in the case of nuclear matter. 



7.1 Infinite Nuclear Matter 

In symmetric, infinite nuclear matter the sources of the fields are constant and can be 
related to the nucleon Fermi momentum kp according to || 



4 r k F ^ 2k 3 
mH{¥)\A) = —J dk = ^ (124) 



4 r k F - M* N 



{A\W(r)\A) = j—- 3 / dk JV , , (125) 

( 27r ) J Jm* n 2 + k 2 



where denotes the constant value of the effective nucleon mass defined by Eq . ([43D , or 
equivalently, Eq.(116). Obviously these equations can also be deduced from the definitions 
(Eqs. flHHD and ( |102| )) of the scalar and vector densities provided that, for p s , one uses 
Eqs.(||) and @ to write 

M " - 1 - M - (126) 



E-V coshe X / M tf + p 



Let (<j,uj) be the constant mean- values of the meson fields. From Eqs.( |109| ) and (|110|) 
we find 

u = — — , (127) 

m 2 „ 
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Figure 1: Mean- field values of the a meson for various bag radii as a function of pb- The 
solid, dotted and dashed curves show g a a for R° B = 0.6, 0.8 and 1.0 fm, respectively. The 
quark mass is chosen to be 5 MeV. 



a 



9a 



m 



C(a) 



(27T) 



dk 



M 



N 



Mtf + k 2 



(128) 



where C(a) is now the constant value of C in the scalar fleld.B As emphasised by Saito 
and Thomas |16fl , the self-consistency equation for a, Eq.( |128|) , is the same as that in 
QHD except that in the latter model one has C(w) = 1 (i.e. the quark mass is infinitely 
heavy) . 

Once the self-consistency equation for a has been solved, one can evaluate the energy 
per nucleon. From the Dirac equation (|119 ), or simply using Eq . (|58|) , the energy of a 
nucleon with momentum k is 



E(k) = V + M* N cosh£ = gju + \J M* N 2 + k 



(129) 



2 Note the change in notation from the earlier papers of Saito and Thomas where C was used for what 
we now call S. 
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Figure 2: Effective nucleoli mass (m 9 =5MeV). The curves are labelled as in Fig.l. 



which contributes to the energy per nucleon by the amount 

4 r k F 



E nucl./ A 



Pb(2tt) 3 

1 



PB 



dkE(k), 

4 fk F 



(27T) 



dk\/Mtf + k 2 



(130) 
(131) 



The contribution of the energy stored in the meson fields is, from Eq.( |113| ), 

1 



E r 



2 Pl 



/ 2—2 2 — 2\ 



(132) 



and using the expression for uj we finally obtain the following expression for the energy 
per nucleon 



j^total I ^ 



PB 



(2k) 



(133) 



We determine the coupling constants, g a and g u , so as to fit the binding energy (—15.7 
MeV) per nucleon and the saturation density (po = 0.15 fm -3 ) for symmetric nuclear 
matter at equilibrium. The coupling constants and some calculated properties of nuclear 



28 



o 



0.9 - 



0.8 - 



0.7 



0.6 



0.5 




I... 

50 100 150 200 250 300 350 400 



g o a (MeV) 

Figure 3: Scalar density ratio, C(cf), as a function of g a a {m q = 5 MeV). The curves are 
labelled as in Fig.l. 

matter (with m q = 5 MeV, m a = 550 MeV and m w = 783 MeV) at the saturation density 
are listed in Table |2|. The most notable fact is that the calculated incompressibility, K, is 
well within the experimental range: K pa 200 — 300 MeV [27 . Also our effective nucleon 
mass is much larger than in the case of QHD. In the last two columns of Table |2] we show 
the relative modifications (with respect to their values at zero density) of the bag radius 
and the lowest eigenvalue, x, at saturation density. The changes are not large. In order 
to show the relative insensitivity to the quark mass (as long as it is small) we note that 



M 



N 



(756, 753) MeV at saturation density and K = (278, 281) MeV for R° B = 0.8 fm 



and m q = (0, 10) MeV, respectively. 

In Figs.|l| and ^ we show the mean-field values of the a meson and the effective nucleon 
mass in medium, respectively. In both cases their dependence on the bag radius is rather 
weak. The scalar density ratio, C(a), and the ratio of the integral 1(a) to Iq are plotted 
as a function of g a a in FigsS and [§, respectively As a increases, the scalar density ratio 
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Figure 4: The ratio of 1(a) to Iq (m q = 5 MeV). The curves are labelled as in Fig.l. 

decreases linearly, while the ratio, I/Io, gradually increases. Here we note that S(0) = 
(0.4819, 0.4827, 0.4834) and J = (0.2421, 0.3226, 0.4028) fm for R° B = (0.6, 0.8, 1.0) fm, 
respectively. 

It would be very useful to give a simple parametrization for C and I/Iq because they 
are completely controlled by only the strength of the local a field. We can easily see that 
C is well approximated by the linear form: 

C(a) = l-ax (g a a), (134) 

with g a a in MeV and a = (6.6, 8.8, 11) xlO" 4 , for R% = (0.6, 0.8, 1.0) fm, respectively. 
For I/Iq we find a quadratic form: 

= 1 + 6i x (g a a) -b 2 x (g a a) 2 , (135) 

Jo 

with bx = (3.7, 4.9, 6.1) xlO" 4 and b 2 = (3.9, 5.2, 6.5) xl0~ 7 . More comments and 
discussion of the results for nuclear matter can be found in the previous publications [|l|, 
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Table 2: Coupling constants and calculated nucleon properties in symmetric nuclear mat- 
ter at normal nuclear matter density. The effective nucleon mass, M^, and the nuclear 
incompressibility, K, are quoted in MeV. The bottom row is for QHD. 



R%(im) 




S2/47T 


M* N 


K 


SR B 


Sx 
XO 


0.6 


5.86 


6.34 


729 


295 


-0.02 


-0.13 


0.8 


5.40 


5.31 


754 


280 


-0.02 


-0.16 


1.0 


5.07 


4.56 


773 


267 


-0.02 


-0.21 


QHD 


7.29 


10.8 


522 


540 







As a practical matter, we note that Eq. (|123|) is easily solved for g a (o~) in the case 
where C(a) is linear in g a a - as we found in Eq. (|134|) . In fact, it is easy to show that 



M* N = M 



N 



(136) 



(recall g a = g a (a = 0), Eq. ( |112j )) so that the effective cxN coupling constant decreases at 
half the rate of C(a). (Equation (|136|) is quite accurate up to twice nuclear matter density.) 
Having explicitly solved the nuclear matter problem by self-consistently solving for the 
quark wave functions in the bag in the mean scalar field one can solve for the properties 
of finite nuclei without explicit reference to the internal structure of the nucleon. All one 
needs is Eqs.( |134[ ) and ( |136| ) for C(a) and as a function of g a a. 
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7.2 Initial results for finite nuclei 

In order to keep the present work to a reasonable length we propose to present detailed 
numerical studies of the properties of finite nuclei calculated within the present framework 
in a later paper . However, to illustrate that the approach has some promise, we have 
performed some preliminary calculations for the doubly closed shell nucleus 16 O. This 
requires the self-consistent solution of our equations for 6 Dirac orbitals - the lsi/ 2 , lpi/2 
and 1^3/2 states of the protons and neutrons. The central Coulomb interaction has been 
taken into account for the protons. 

The numerical calculation was carried out using the techniques described by Walecka 
and Serot ||. The resulting charge density for 16 is shown in Fig.|5] (dotted curve) 
in comparison with the experimental data (hatched area) and QHD ||. In this 

o. 10 
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^ 0.06 

I 

o." 0. 04 

0.02 
0.00 

1 2 3 4 5 

r (fm) 

Figure 5: The charge density of w O in the present model and QHD, 
compared with the experimental distribution. 

calculation we used the parameters given in Table |2| for R° B = 0.8 fm. However, as the 
central density of 16 was a little high we increased the model-dependent slope of the scalar 
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Figure 6: Scalar, vector and coulomb potentials for 16 in the present 
model. 

density C(a) and the coupling constant g a (cr) (i.e. the parameter a in Eqs.( |134j ) and (|136|) 
) by 10% to obtain the result shown. The corresponding effect on the saturation energy 
and density of nuclear matter is very small: po — > 0.1496 fm~ 3 and the energy per nucleon 
becomes —15.65 MeV. In Fig.^, we also show the scalar and vector fields corresponding 
to the charge density of 16 that was shown in Fig.[5[ Using the same parameter set one 
also finds a very reasonable fit to the charge density of 40 Ca. It is presented in Fig.[7]. 

To conclude this initial investigation of finite nuclei we record, in Table ^|, the calcu- 
lated single particle binding energies of the protons and neutrons in A0 Ca, in comparison 
with the results of QHD and the experimental data ||. Because of the smaller scalar and 
vector field strengths in the present model, compared with QHD, the spin-orbit splitting 
tends also to be smaller - perhaps only | of the experimentally observed splittings. In this 
context it is interesting to show the sensitivity of the present model to just one feature 
of the underlying structure of the nucleon, namely the mass of the confined quark. As 
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Figure 7: The charge density of 40 Ca in the present model and QHD, 
compared with the experimental distribution. 

an example, we consider the case m q = 300 MeV, which is a typical constituent quark 
mass. For m q = 300 MeV and Ro = 0.8 fm, the coupling constants required to fit the 
saturation properties of nuclear matter are <j£/47r = 6.84 and g^,/47r = 8.51, and the ef- 
fective nucleon mass (at saturation) and the incompressibility become 674 MeV and 334 
MeV, respectively. Using these parameters one finds that the charge densities of 16 and 
40 Ca are again reproduced very well (without any need to vary the slope parameter "a", 
i.e., a = 3.9 x 10 -4 ). From the table, one can see that a heavier quark mass gives a 
spectrum closer to that of QHD - as discussed by Saito and Thomas in Ref. [FJ| - and 
in better agreement with the observed spin-orbit splittings. In conclusion we notice one 
more point: the 2s\/2 and lc/3/2 levels are inverted compared with the experimental data 
(see also Ref.||). This may be connected with the effect of rearrangement [III 50), which 
is not considered here. Also, the effect of antisymmetrization remains to be investigated. 
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Table 3: Predicted proton and neutron spectra of 40 C a compared with QHD and the 
experimental data. QMC(5, 300) means the present model with Rq = 0.8 fm and m q = 5 
and 300 MeV, respectively (see text). All energies are in MeV. 



Shell 




neutron 






] 


proton 




QHD 


QMC(5) 


QMC(300) 


Expt. 


QHD 


QMC(5) 


QMC(300) 


Expt. 


lSl/2 


54.9 


43.5 


47.1 


51.9 


46.7 


35.5 


38.9 


50±10 


1P3/2 


38.6 


32.5 


34.5 


36.6 


30.8 


24.7 


26.7 


34±6 


lPl/2 


33.1 


31.3 


32.4 


34.5 


25.3 


23.5 


24.5 


34±6 


14/2 


22.6 


20.0 


21.0 


21.6 


15.2 


12.6 


13.6 


15.5 


2si/ 2 


14.6 


15.0 


15.3 


18.9 


7.4 


7.7 


7.9 


10.9 


14/2 


14.1 


17.9 


17.3 


18.4 


6.8 


10.5 


9.8 


8.3 
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8 Summary 



Starting with a hybrid model in which quarks confined in nucleon bags interact through 
the exchange of scalar and vector mesons, we have shown that the Born-Oppenheimer 
approximation leads naturally to a generalisation of QHD with a density dependent scalar 
coupling. The physical origin of this density dependence, which provides a new saturation 
mechanism for nuclear matter, is the relatively rapid increase of the lower Dirac component 
of the wavefunction of the confined, light quark. We confirm the original discovery of 



Guichon |L5[ that once the scalar and vector coupling constants are chosen to fit the 
observed saturation properties of nuclear matter the extra, internal degrees of freedom 
lead to an incompressibility that is consistent with experiment. 

In the case of finite nuclei we have derived a set of coupled differential equations 
which must be solved self-consistently but which are not much more difficult to solve 
than the relativistic Hartree equations of QHD. Initial results for 16 (see Fig.|5|) are quite 
promising but a full numerical study will be presented in a later work [p8fl . 

The successful generalisation of the quark-meson coupling model to finite nuclei opens 
a tremendous number of opportunities for further work. For example, earlier results for the 
Okamoto-Nolen-Schiffer anomaly [20|, the nuclear EMC effect |[22|| , the charge-symmetry 
violating correction to super-allowed Fermi beta-decay |2l| and so on, can now be treated 
in a truly quantitative way. 

It will be very interesting to explore the connection between the density dependence of 
the variation of the effective cx-nucleon coupling constant, which arises so naturally here, 
with the variation found empirically in earlier work. We note, in particular, that while 
our numerical results depend on the particular model chosen here (namely, the MIT bag 
model), the qualitative features which we find (such as the density dependent decrease of 
the scalar coupling) will apply in any model in which the nucleon contains light quarks and 
the attractive N-N force is a Lorentz scalar. Of course, it will be important to investigate 
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the degree of variation in the numerical results for other models of nucleon structure. 

We could list many other directions for future theoretical work: for the replacement 
of the MIT bag by a model respecting PCAC (e.g. the cloudy bag model |nj]), the re- 
placement of ex-exchange by two-pion exchange, the replacement of u exchange by nucleon 
overlap at short distance, the inclusion of the density dependence of the meson masses 
themselves 117 , |32| and so on. On the practical side, we stress that the present model can 
be applied to all the problems for which QHD has proven so attractive, with very little 
extra effort. It will also be interesting to explore its phenomenological consequences in 
this way. 
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Appendix 

Here we want to demonstrate that the center of mass correction to the bag energy is 
essentially independent of the external scalar field. (Since the vector fields do not alter 
the quark structure of the nucleon we need not consider them.) To avoid the difficulties 
associated with the confinement by a sharp boundary we consider a model where the 
quark mass grows quadratically with the distance from the center of the bag. This is 
justified because we do not look for the cm. correction itself but only for its dependence 
on the external field. Moreover, after the strength of the confining mass has been adjusted 
to reproduce the lowest eigenfrequency of the bag, we have found that the corresponding 
wave functions are rather similar to those for the bag. 

To estimate the cm. correction we also make the assumption that the quark number is 
a good quantum number, which allows us to formulate the problem in the first quantized 
form. This amounts to neglecting the effect of quark- ant iquark excitation and is therefore 
not a very strong constraint. 

Thus the model is defined by the following first-quantized Hamiltonian: 

Hb = Yl 7o(«)[t(*) -Pi + m(fi)], Pi = -iVi, (137) 

i=l,N 

with 

m(f) =m* + Kr 2 , (138) 

where m* is the mass of the quark in the presence of the external scalar field. 

By assumption, N is a number, so we can define intrinsic coordinates (p, n) by 

— * 

P -> - 

= Pi~Jj' P = HPii J2^i = ^ ( 139 ) 

i i 

p t = r t -R, i?=l^rl, 5>< = °- (140) 

i i 

Then we can write the Hamiltonian in the form: 

Hb = H-intr. + H C m, (141) 
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H intr . = E7o(i)[7(i)-^ + "i(Pi)], (142) 

i 

— * 

H cm = ^ ■ J^ToWtW + J2^(i)[m(fi) - m{pi)\. (143) 

i i 

This separation into an intrinsic and a cm. Hamiltonian is correct because: 

— * — * 

1. H intr commutes with P and R, 

2. One has [H C m,R] = — i Hi 7o(*)7(0- Since for a Dirac particle 707 is the velocity, 
one can identify the RHS of the previous equation with the time derivative of R, 
which is consistent. 

The fact that the cm. Hamiltonian depends on the intrinsic coordinates is not a surprise 
because the separation is only complete in certain special cases. 
We now look for the intrinsic energy of the bag, writing 

Hintr. = H B — H C m- (144) 

All that we know are the eigenstates of Hb but we can consider Hcm as a correction of 
order 1/JV with respect to the leading term in the bag energy. Therefore we estimate its 
effect in first order perturbation theory, that is 

Hintr. = E B — (B\H C m\B) = E B — E C m, (145) 
where \B) is the eigenstate of H B with energy E B . We must therefore evaluate 

— * 

Ecu = (B\^ ■ 5>oO')70') + 5>o(i) [m(rl) - m(f t - R)] \B) 

j i 

= ■ £7oO')7(i) + 2KR .£70(0* - KR^^IB). (146) 

j i i 

(Note that there are no ordering problems as long as iV is a number). Let \ct) be the one 
body solutions, that is (in units such that R B — 1) 

7o[7"P + m(r)}(f) a = Vt a <p a . (147) 
39 



If we assume that \B) has all the quarks in the lowest mode then elementary techniques 
for many-body systems lead to the results 

(B\^ -J2^Um)\B) = fi -(0|7oK + ^r 2 )|0), 
j 

(B\R-J2lom\B) = (0| 7 or 2 |0), 

i 

(B\R"J2lo^\B) = ^(0| 7 or 2 |0)+(l-l)(0| 7o |0)(0|r 2 |0). 



(148) 



so that we get 

E CM = fi -m*(0| 7 o|0)+^(l-^) ((0| 7o r 2 |0)-(0|7o|0)(0|r 2 |0) 



N 

Q - m*(0| 7 o|0> + K ((0| 7o r 2 |0) - (0| 7o |0) (0|r 2 |0)) + 0(1/N) 



(149) 



Note that we keep only the leading term in 1/N. This is consistent with our initial 
approximation according to which Eqm is computed as a correction of order 1/N with 
respect to the leading term in the bag energy. 

To proceed we need to evaluate the single particle matrix elements which appear in 
the expression for Ecu- To determine the wave function we solve Eq.( |147] ) numerically 
and adjust the constant K to give f2 = 2.04 - i.e. the lowest energy level of the free bag, 
in units such that Rb = 1- We found K = 1.74. 

Then we compute Eqm numerically according to Eq. ( |149|) , as a function of m*. The 



result is shown in Fig|| where we also plot the value of Qq- One can see that in the range 
— 1.5 < m* < 0, which certainly contains the possible values of m* in the case of finite real 
nuclei, the value of Eqm is almost constant. For instance at m* = —1, Ecu differs from 
its free value by only 6%. Furthermore, its variation is clearly negligible with respect to 
that of Qq. Thus, for practical purposes, it is a very reasonable approximation to ignore 
the dependence of Eqm on the external field. 
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Figure 8: Dependence of Eqm (full line) and Qq (dashed line) on m* 
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